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ABSTRACT 

Anticipating forthcoming Chandra X-ray observations, we describe the continuation of inter- 
stellar cooling flows deep into the cores of elliptical galaxies. Interstellar gas heated to T > 1 
keV in the potential of massive black holes (r < 50 pc) should be visible unless thermal heating 
is diluted by non-thermal pressure. Since our flows are subsonic near the massive black holes, 
distributed cooling continues within '^ 300 pc. Dark, low mass stars formed in this region may 
be responsible for some of the mass attributed to central black holes. 

Subject headings: galaxies: elliptical and lenticular - galaxies: cooling flows - galaxies: interstellar 
medium - galaxies: black holes 



1. INTRODUCTION 

Considerable observational evidence now exists 
for supermassive black holes (SMBHs) in the cen- 
ters of giant elliptical galaxies. These mass con- 
centrations, long suspected to be the origin of non- 
thermal energy in active galaxies, are now thought 
to be responsible for large gas or stellar velocities 
in galactic cores (Ho 1998). We discuss here the 
interaction of SMBHs with hot (T - lO"^ K) in- 
terstellar gas. The local compression and heating 
of interstellar gas in the SMBH potential can be 
detected with the Chandra X-ray telescope. In- 
deed, the inflow of interstellar gas may explain the 
growth of SMBHs in galactic cores. However, X- 
ray detection of SMBHs may not be possible if the 
non-thermal interstellar pressure is comparable to 
the pressure of hot thermal gas. 

As we were completing this work, a similar 
study of the SMBH-gas interaction by Quataert 
& Narayan (1999) (QN) appeared on the astro-ph 
preprint server. While some of the conclusions of 
QN are supported by our flow calculations, our re- 
sults are qualitatively different in a number of im- 



portant ways. The steady state flows considered 
by QN are constrained to pass through a sonic 
point, becoming supersonic near the SMBH. By 
comparison, our (non-steady) computational gas 
flows near the SMBH adapt naturally to the in- 
terstellar cooling flow in the surrounding galaxy 
and, as a result, our flow velocity is fully sub- 
sonic near the SMBH. We also consider several 
additional physical effects not discussed by QN: 
time-dependent ejection of mass by evolving stars, 
heating by Type la supernovae, non-thermal pres- 
sure, and the gravitational field of (dark) cooled 
interstellar gas which in some observations may be 
confused with the mass of the SMBH. 

Although we use observations of the bright 
Virgo elliptical NGC 4472 as a guide for our mod- 
els, we do not seek here a set of parameters for an 
optimal fit to this galaxy. Instead, our objective is 
to illustrate some of the physical effects that must 
be considered in interpreting X-ray observations 
of the cores of bright ellipticals with the Chandra 
telescope. At a distance d = 17 Mpc appropriate 
for NGC 4472, Chandra's resolution (^ 0.5 arcsec- 
onds) is ~ 40 pc. 



2. CONSTRUCTING MODELS 

The computational procedure wc follow here is 
identical to that described in our recent papers 
(e.g. Brighenti & Mathews 1999a). With NGC 
4472 in mind, we construct a galaxy in a cosmo- 
logical environment beginning with an overden- 
sity perturbation designed to form a galaxy group. 
Stars of total mass M^t = 7.26 x 10" Mq (appro- 
priate to NGC 4472 with M/Lb = 9.2; van der 
Marel 1991), are formed from the baryons at an 
early time (t^s ~ 1 Gyr) and later (t* = 2 Gyrs) 
merge into a de Vaucouleurs profile. An evolving 
NFW dark halo develops around the galaxy over 
time, accompanied by secondary infall, although 
within an effective radius (rg — 8.57 kpc for NGC 
4472) stellar mass loss provides about half of the 
hot interstellar gas. A cooling flow develops in 
the hot interstellar gas at early times and slowly 
evolves toward its present configuration at t„ = 13 
Gyrs. Over time a large mass of interstellar gas 
cools in NGC 4472, Mdo « tnL^/i5kT/2nmp) « 
4 X 10^° Mq. Since this mass exceeds the mass 
observed in the central black hole in NGC 4472, 
~ 2.6 X 10^ Mq (Magorrian et al. 1998), inter- 
stellar gas must cool in localized thermal insta- 
bilities distributed throughout an extended galac- 
tic volume within Tg. This distributed interstellar 
cooling is modeled by adding a sink term to the 
equation of continuity: —q{r)p/tdo oc p^, where 
tdo = 5mpkT/2iipA is the time for gas to cool lo- 
cally by radiative losses at constant pressure (e.g. 
Sarazin & Ashe 1989). We find that taking q = 1 
provides a satisfactory global interstellar flow so- 
lution for NGC 4472 that matches currently ob- 
served density, temperature and metallicity pro- 
files reasonably well (Brighenti & Mathews 1999a; 
1999b); we adopt q = 1 here. 

To resolve the small scale flow near the SMBH 
we have refined our logarithmic spatial grid in 
the galactic core to an innermost zone of ra- 
dius 12 pc. Interstellar gas is influenced by the 
black hole within a radius rh ^ GMth/c^ ~ 
48(Af/lO^M0)(c2/3OO km/s)-^ pc; c^ is the sound 
speed in the hot gas. We ignore the effects of 
galactic rotation. Although the mass-losing stel- 
lar systems in large ellipticals like NGC 4472 are 
slowly rotating, their X-ray images do not indicate 
the flattening expected if the hot interstellar gas 
shared this rotation (Brighenti & Mathews 1996). 



Evidently, efficient cooling dropout or turbulent 
viscosity removes angular momentum from the hot 
gas. 

3. A "STANDARD" MODEL 

Figure 1 shows the interstellar density, temper- 
ature and (ROSAT band) X-ray surface brightness 
evolved to the present time, i„ — 13 Gyrs. Dur- 
ing the evolution, new gas and energy are supplied 
from stellar sources: (1) the specific stellar mass 
loss rate is a^{t) = 4.7 x 10-'^°[t/{ta - t,^)]"^-^ 
sec~^ and (2) gas is heated by Type la super- 
novae at a rate SNu(i) :=SNu(t„)(i„/i), where 
SNu(t„) = 0.03 SNIa per 100 yrs per 10^° Lbq 
(Cappellaro et al. 1997). The interstellar gas tem- 
perature is nearly virial, T « 10^ K. The stellar 
density has a de Vaucouleurs profile but flattens 
to /9, oc r^"'® within a "break" radius ri, ~ 206 pc 
(Faber et al. 1997). 

Dashed lines in Figure 1 refer to solutions with- 
out central black hole and the solid lines refer to 
solutions with a central black hole of mass 2.0 x 10^ 
Mq introduced at t, = 2 Gyrs. In this q ~ I so- 
lution, an additional 0.7 x 10® Mq flows inside 
r = 12 pc by t„ = 13 Gyrs, increasing the cen- 
tral dark mass to 2.7 x 10" Mq. The gas density 
(with SMBH) varies approximately as n ex r~^'^ 
over most of the flow. For r < 300 pc, the com- 
puted density appears to exceed the innermost ob- 
servations; we discuss this apparent discrepancy 
in more detail below. Light curves in Figure 1 for 
T{r) refer to the background flow; heavy curves 
show the projected temperature including emis- 
sion from distributed cooling regions and corre- 
spond to the temperature observed. For constant 
q, cooling dropout is important throughout the 
flow and extends to the region adjacent to the 
SMBH. The gas temperature rises dramatically 
within about 50 pc from the SMBH - detection 
of this thermal peak by Chandra would provide 
additional evidence for a SMBH. This tempera- 
ture spike should be visible even when viewed with 
Chandra in projection against the entire cooling 
flow: the mean projected thermal temperatures 
within apertures oi R = 50, 100 and 500 pc are 
2.95, 2.70 and 2.00 xlO^ K (with SMBH) and 
1.04, 1.05 and 1.09 xlO^ K (without SMBH) re- 
spectively. 

The upper panel in Figure 2 shows the cur- 



rent flow velocity v(r) and the local mass flow 
M{r) = Airr^pv. A small influence due to the 
SMBH potential at times t ^ t,: can be seen in 
the distant flow r > 30 kpc where the flow time 
tf — "''1'^ ~ in — i*- Overall, the flow is very 
subsonic even in the vicinity of the SMBH where 
Ca ^ 300 km s^^. This occurs for two reasons: (1) 
the high pressure in the central thermal spike (e.g. 
Holzer & Axford 1970), enhanced by the SMBH, 
and (2) mass dropout which tends to make sub- 
sonic flow even slower. By contrast, the solutions 
of QN (with a*(t) = 0) are constrained to pass 
through a sonic transition, causing the gas to ac- 
celerate through the thermal spike. [Without mass 
dropout (g = 0) we also find a transition to super- 
sonic flow at ^ 60 pc] Our calculations shown in 
Figures 1 and 2 were done with boundary condi- 
tion u{r = 0) = 0; accumulated gas in the cen- 
tral zone cools there. To encourage the onset of 
supersonic flow near the SMBH, we repeated the 
calculation allowing the gas to flow freely through 
a small radius (as small as r = 1 pc with a very 
high resolution grid) and found that the flow was 
essentially unchanged from that shown in the Fig- 
ures. The fully subsonic character of our solu- 
tions is therefore quite robust and independent of 
our boundary condition at r = 0. [In principle 
our flow should become supersonic very near the 
Schwarzschild radius (jg w 10""* pc), but this can- 
not be resolved by our spatial grid and it is likely 
that additional physical effects become important 
as T approaches r^.] 

Because our flow in the thermal spike near the 
SMBH is subsonic, the cooling time t^o remains 
comparable to the flow time t/;o«) — t jv through- 
out this region and cooling dropout continues un- 
abated at small radii, unlike the QN solutions. 
This concentrated mass dropout has two impor- 
tant astrophysical consequences: (1) the mass de- 
posited by interstellar cooling near the SMBH may 
be comparable to that of the SMBH itself, and (2) 
the rate that mass flows into the SMBH at the ori- 
gin is significantly reduced (by about ^ 10) below 
the QN solutions. 

The mass of cooled interstellar gas in our model 
with SMBH in Figure 2 is large: Ms.o{r) = 2.3 x 
10^, 3.5 X 10^ and 8.0 x 10^ Mq within r = 50, 
100 and 300 parsecs respectively. These masses 
are comparable to SMBH masses found by Magor- 
rian et al. (1998) and they occur within the 



small region (^ 4 arcseconds or 330 pc for NGC 
4472) where stellar velocities are enhanced by the 
SMBH. The gravitational acceleration due to the 
SMBH, the dropout mass, old stars and dark halo 
are also shown in Figure 2; note that the dropout 
mass dominates gravitational forces at r ~ 200 
pc. In view of the possibility of contamination 
with dropout mass, the true mass of the central 
SMBH may be lower than those quoted by Magor- 
rian et al. These considerations depend on our as- 
sumption that interstellar mass dropout continues 
to the smallest galactic radii; at the present time 
there is no reason to question this assumption. 

In Figure 3 we illustrate a solution identical to 
that in Figure 2 but with an initial SMBH mass 
of only M\,h = 1 x 10^ M©; the visibility of the 
thermal spike is noticeably reduced. 

4. "NON-STANDARD" MODELS 

We now return to the apparent disagreement in 
Figure 1 between our models and the innermost 
observed Sa;(i?) and deprojected electron density. 

Can the apparent flattening of observed S^, (i?) 
at small radii be due to absorption of X-rays by 
cold gas in r < 300 pc? The HII gas that pro- 
duces optical emission lines in NGC 4472 (Mac- 
chetto et al. 1996) cannot account for the X-ray 
absorption since H/3 emission from such HII clouds 
in pressure equilibrium would be > 1000 times 
more intense than L/f/j observed in this region of 
NGC 4472. Suppose instead that the absorbing 
gas is cold HI or H2 and is distributed in a uni- 
form disk extending to r ^ 300 pc oriented per- 
pendicular to the line of sight. If the disk column 
density is iV w 10^^ cm'^ (i.e. r(l keV) « 1), 
the total disk mass, ^ 2 x 10^ M©, would be 
consistent with upper limits on observed HI and 
H2 in NGC 4472, M^oU ^ 10^ Mq (Bregman, 
Roberts & Giovanelh 1988; Braine, Henkcl & Wik- 
lind 1997). However, in pressure equilibrium with 
the hot cooling flow gas, the density of such cold 
gas (at T - 15 K Mathews & Brighenti 1999) is 
very high, Jicou ^ 6x10^ cm~^, and the disk thick- 
ness is incredibly small, h « N/ricoid ~^ 6 x 10^^ 
pc. Whether such a peculiar and fragile disk 
could exist in NGC 4472 (and other Virgo ellip- 
ticals with similar cores; Brighenti & Mathews 
1997) seems unlikely. However, optical line emis- 
sion from the surface of the disk ionized by stel- 



lar UV (Mathews & Brighenti 1999) would exceed 
observed limits. For HII gas at T « 10^ K the 
pressure equilibrium density at rd — 300 pc is 
uhii ~ 850 cm~'^ and the Stromgren radius is 
rs ~ 10^® cm. The total H/3 energy emitted by 
the disk, LHfi.d « 2TTrlrsn'jjjjeHi3 « 4 x 10^^ ergs 
s~^, is '-~^ 6 times larger than the total H/3 luminos- 
ity from NGC 4472, LHf3 « 7 x lO^^ ergs s'^. This 
inconsistency cannot be removed with other geo- 
metrical arrangements of the absorbing cold gas 
within 300 pc, so the apparent inner flattening in 
Figure 1 is probably not due to absorption. 

We have also explored the possibility that the 
T,x{R) observations shown in Figure 1 were cen- 
trally flattened by the 4 arcsecond half power 
Gaussian PRF (Giaconni et al. 1979) of the Ein- 
stein HRI. However, such a convolution of our 
computed T^xiR) cannot duplicate the core ob- 
served in NGC 4472. X-ray cores in NGC 4472, 
NGC 4636 and NGC 4649 (Trincheri, Fabbiano & 
Canizares 1986) are not artifacts of limited obser- 
vational resolution. 

These difficulties led us to explore two ad- 
ditional "non-standard" models for NGC 4472: 
flows (1) with extreme mass dropout inside 300 
pc or (2) with additional non-thermal pressure in 
this same region. 

It is often claimed that mass dropout reduces 
the X-ray surface brightness in the central regions 
of cooling flows, but this is not true in general. 
For flows in which most of the gas in the central 
region has advected inward from large radii, cool- 
ing dropout at large radii can reduce the amount 
of gas and associated emission T,x{R) in the cen- 
tral core. However, this does not apply to galac- 
tic cooling flows in which much of the gas in the 
central flow originates from stars within or near 
this region. In addition, the X-ray emissivity and 
T,x{R) are locally enhanced at (denser) coohng 
sites where mass dropout occurs. To illustrate the 
futility of lowering the central T,^{R) in NGC 4472 
with increased local cooling, we show in Figure 3 
density and temperature profiles for the standard 
flow (at tn and with SMBH) but with a cooling 
dropout that truncates the background flow den- 
sity: q{r) = e"('')/"'" with n„ = 0.006 cm^^. The 
thermal spike around the SMBH is still present in 
this solution, but the apparent n{r) (and T,x{R)) 
within several 100 pc is only slightly reduced. 
The mean projected temperature within 100 pc 



is T{R < 100) == 1.9 X 10^ K. 

Finally, we consider cooling flow models in 
which n{r) and T.x{R) in the inner flow arc lowered 
by an additional non-thermal pressure due to rel- 
ativistic particles or magnetic flclds, as suggested 
by radio emission observed in most elliptical cores. 
To explore this possibility we show in Figure 3 the 
standard flow (at i„ = 13 Gyrs with SMBH) in 
which an additional pressure is introduced in the 
galactic core: Ptot = -P(l + Se^*"/^ '^p^), where P 
is the thermal gas pressure. The density (and T.^) 
in the cooling flow within ~ 300 pc are reduced by 
this additional pressure and parameters for Ptot (r) 
could be adjusted in an ad hoc fashion to achieve 
a more perfect fit with the observations. However, 
the apparent temperature of the thermal gas is 
dramatically lowered since the external flow is now 
supported mostly by non-thermal pressure. The 
projected apparent temperature within R = 100 
pc is only T = 0.40 x 10^ K, assuming this gas 
is not heated by collective interactions with the 
relativistic gas. The discernibility of the thermal 
spike near the SMBH is now greatly reduced. 

5. FINAL REMARKS AND CONCLU- 
SIONS 

In this study we have examined central cool- 
ing flows in a luminous elliptical similar to NGC 
4472. We have emphasized the evolution of qui- 
escent flows that arise naturally from gas and en- 
ergy supplied by galactic stars and by secondary 
infall in the outer regions. We have not explored 
possible transient interactions between the cool- 
ing flow and an intermittently active galactic core 
(e.g. Ciotti & Ostriker 1997). Our objective here 
has been to anticipate some of the interpretive is- 
sues that may arise when Chandra observations 
become available. 

Our main conclusions arc as follows: 
(1) Some dark matter currently attributed to the 
central supermassive black hole (SMBH) may be 
due to a population of non- luminous stars created 
from cooled interstellar gas having a bottom-heavy 
IMF. In the standard coohng flow for NGC 4472 
with SMBH, neutral clouds become gravitation- 
ally unstable (and truncate the IMF) at masses 
ruu = 0.2, 0.3 and 0.8 Af© at r = 50, 100 and 
300 pc respectively (Mathews & Brighenti 1999). 
Such dark stars, possibly having approximately ra- 



dial orbits, can also produce a radial variation in 
the global baryonic M /L. 

(2) A successful observation of a high tempera- 
ture peak in r < 50 pc by Chandra in a nearby 
bright elliptical would provide additional strong 
evidence for the existence of SMBHs. However, if 
additional non-thermal pressure is present in ellip- 
tical cores the SMBH thermal peak may no longer 
be discernible. 

(3) Our central flows (r < 300 pc) differ from the 
Bondi-type transonic flow described by Quataert 
& Narayan (1999) since we do not insist on a sonic 
transition. The mass flow rate into the SMBH is 
considerably reduced in our models. 

(4) As long as our standard cooling flow assump- 
tions hold, we expect the gas density to rise 
steadily toward the SMBH approximately as n oc 
j,-5/4_ g^^ current observations suggest that such 
a gas density peak is not present. About half of 
the gas within 300 pc in NGC 4472 is supplied by 
mass loss from stars in this region. In our standard 
cooling flow model we have assumed that all stel- 
lar ejecta rapidly enters the hot interstellar phase 
(in < 10^ yrs). It has occasionally been suggested 
(e.g. Thomas 1986) that steUar ejecta cools be- 
fore entering the hot gas, effectively reducing a*. 
However, we find that the central peak in n{r) and 
Sa; {R) persists even if the stellar mass loss rate a* 
is reduced at small r. Moreover, dynamical argu- 
ments and gas temperature and metallicity gra- 
dients observed at larger radii (r < 3re) strongly 
indicate that cooler, enriched stellar ejecta is con- 
ductively melting into the hot phase. 

Studies of the evolution of hot gas in ellipti- 
cal galaxies at UC Santa Cruz are supported by 
NASA grant NAG 5-3060 and NSF grant AST- 
9802994 for which we are very grateful. FB is 
supported in part by Grant MURST-Cofin 98. 

REFERENCES 

Braine, J., Henkel, C. & Wikhnd, T. 1997, A&A 
321, 765 

Bregman, J. N., Roberts, M. S. & Giovanelli, R. 
1988, ApJ 330, L93 

Brighenti, F. & Mathews, W. G. 1996, ApJ 470, 
747 



Brighenti, F. & Mathews, W. G. 1997, ApJ 486, 
L83 

Brighenti, F. & Mathews, W. G. 1999a, ApJ 512. 
65 

Brighenti, F. & Mathews, W. G. 1998b, ApJ (in 
press) 

Cappellaro, E., Turatto, M, Tsvetkov, D. Yu., 
Bartunov, O. S., PoUas, C., Evans, R. & 
Hamuy, M. 1997, 1997, Astron & Ap. 322, 
431 

Ciotti, L. & Ostriker, J. P. 1997, ApJ 487, L105 

Faber, S. M. et al. 1997, AJ, 114, 1771 

Giacconi, R. et al. 1979, ApJ, 230, 540 

Ho, L. C. 1999, in Observational Evidence for 
Black Holes in the Universe^ ed. S. K. 
Chakrabarti (Dordrecht: Kluwer), p 000 

Holzer, T. E. & Axford, W. I. 1970, ARA&A 8, 
31 

Irwin, J. A. & Sarazin, C. L. 1996, ApJ 471, 683 

Macchctto, F., et al. 1996, A&AS 120, 463 

Magorrian, J., et al. 1998, ApJ 115, 2285 

Mathews, W. G. & Brighenti, F. 1999, ApJ (in 
press) (astro-ph/9907364) 

Quataert, E. & Narayan, R. 1999, ApJ (in press) 
(astro-ph/9908199) 

Sarazin, C. L. & Ashe, G. A. 1989, ApJ, 345, 22 

Thomas, P. A. 1986, MNRAS 220, 949 

Trinchieri, G., Fabbiano, G., & Canizares, C. R. 
1986, ApJ, 310, 637 

van der Marel, R. P. 1991, MNRAS, 253, 710 

Fig. 1. — Structure of the standard cooling flow 
solution, shown both with {solid line) and without 
{dashed line) a supermassive black hole (SMBH) of 
mass Mbh = 2 x 10^ Mq. Top Panel: Radial vari- 
ation of apparent gas density including cooling re- 
gions. Middle Panel: Variation of gas temperature 
with galactic radius in the background flow {light 
lines) and with projected radius including emis- 
sion from cooling regions {heavy lines). Bottom 
Panel: Variation of ROSAT band X-ray surface 
brightness with projected radius. Filled circles are 



This 2-column preprint was prepared with the AAS lAT^TJX 
macros v5.0. 



Einstein HRI observations (Trinchieri, Fabbiano & 
Canizarcs 1986) and open circles are ROSAT HRI 
and PSPC observations (Irwin & Sarazin 1996). 

Fig. 2. — Additional properties of the standard 
cooling flow solution shown in Fig. 1. Top Panel: 
Radial variation of the flow velocity (light lines) 
and mass flow rate (heavy lines). Solutions with 
and without supermassive black hole (SMBH) are 
shown as solid and dashed curves respectively. 
Bottom Panel: Gravitational acceleration due to 
SMBH [long dashed line), old stellar population 
[light dotted line), stars formed from cooled inter- 
stellar gas (heavy dotted line), the dark halo (solid 
line) and the interstellar gas (dot-dashed line). 

Fig. 3. — Three additional cooling flow solutions: 
with lower SMBH of mass Mbh = 10^ M© (dotted 
lines); with severe dropout q(r) — exp(n(r)/0.006) 
(dashed lines); and with extra pressure Ptot — 
P(l + 5exp(— r/2kpc) (solid lines). Top Panel: 
Radial variation of apparent gas density including 
cooling regions; for the flow with severe dropout 
the background flow is also shown (light dashed 
line) . Bottom, Panel: Radial variation of gas tem- 
perature in the background flow (light lines) and 
the combined temperature of the background flow 
and cooling regions viewed in projection (heavy 
lines). Nomenclature for the observations is the 
same as in Fig. 1. 
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